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ABSTRACT 

i?-band intensity measurements along the light curve of Type la supernovae 
discovered by the Supernova Cosmology Project (SCP) are fitted in brightness 
to templates allowing a free parameter the time- axis width factor w = s x (1 + 
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z). The data points are then individually aligned in the time-axis, normalized 
and i^-corrected back to the rest frame, after which the nearly 1300 normalized 
intensity measurements are found to lie on a well-determined common rest-frame 
5-band curve which we call the "composite curve". The same procedure is 
applied to 18 low-redshift Calan/Tololo SNe with z < 0.11; these nearly 300 
Z?-band photometry points are found to lie on the composite curve equally well. 
The SCP search technique produces several measurements before maximum light 
for each supernova. We demonstrate that the linear stretch factor, s, which 
parameterizes the light-curve timescale appears independent of z, and applies 
equally well to the declining and rising parts of the light curve. In fact, the B 
band template that best fits this composite curve fits the individual supernova 
photometry data when stretched by a factor s with \ 2 1 DoF ps 1, thus as well as 
any parameterization can, given the current data sets. The measurement of the 
date of explosion, however, is model dependent and not tightly constrained by 
the current data. We also demonstrate the 1 + z light-curve time-axis broadening 
expected from cosmological expansion. This argues strongly against alternative 
explanations, such as tired light, for the redshift of distant objects. 

Subject headings: supernovae: general - cosmology: observations 

1. Introduction 

In systematic searches over the last 12 years the Supernova Cosmology Project (SCP) 
has discovered and studied over 100 supernovae, most of which have been spectroscopically 
identified as Type la supernovae (SNe la) at redshifts up to z = 1.20. Forty-two of these 
events have so far been used for measurements of the cosmological parameters Qm and fl\ 
(Perlmutter et al. (1999), hereafter "P99"; see also Perlmutter et al. 1995 and 1997, hereafter 
"P95" and "P97"). See also Riess et al. 1998 who used sixteen high redshift SNe in a parallel 
study of the cosmological parameters. An important element of these studies has been the 
recognition of homogeneity within a one parameter family of SN la light curves. 

Several approaches have been used to characterize the SN la light-curve family. Phillips 
(1993) first suggested that the range of SNe la light curves might be grouped into a single 
family of curves parameterized by their initial rate of decline (see also the earlier suggestions 
of Pskovskii (1977, 1984)). Phillips (1993) and Hamuy et al. (1995) further observed that the 
absolute magnitudes of SNe la at maximum light were tightly correlated with this decline rate 
in the 5-band light curve — brighter SNe la have wider, slower declining light curves. Riess, 
Press, & Kirshner (1995, 1996) introduced an alternative characterization of this range of 
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light curves, in which a scaled "correction template" is added or subtracted from a standard 
template, effectively widening or narrowing the shoulders of the light curve and thus varying 
its timescale. In our analyses (P95, P97, P99) we instead parameterized the timescale of 
the SN la light curve by a simple "stretch factor", s, which linearly stretches or contracts 
the time axis of a template light curve around the date of maximum light, and thus affects 
both the rising and declining part of the light curve for each supernova. All three methods 
of characterizing the light curve shape or timescale have been used to calibrate the peak 
magnitudes. 

In this paper, we use the dataset of supernovae studied in P99 to examine the empirical 
behavior of SN la light curves, and show that this single parameter, the stretch factor s, 
can effectively describe almost all of the diverse range of £?-band light curve shapes over the 
peak weeks. The stretch factor method was originally introduced and tested (P95) without 
the benefit of either detailed or early data. We can now test whether it applies to the 
entire curve, rather than just the post-maximum decline and determine the scatter about 
this scaled curve. 

This paper is organized as follows. In Sect. 2 we present the set of supernova and the 
data used in this analysis. In Sect. 3 we use fits to a pre-existing template light curve to 
cross-compare our data from each supernova and to produce a single composite plot of all the 
data together. We use the dataset of the composite plot to test the validity of our light-curve 
model and of the stretch parameter. In Sect. 4, we test the hypothesis that redshift is due 
to cosmological expansion by studying the predicted (1 + z) broadening of the light curves 
and we look for evidence for possible supernova evolution. Finally in Sect. 5, we present a 
general discussion. 

2. Data 

A paper based on our first seven SNe gives the details of our photometric and spectro- 
scopic measurements, as well as the analysis methods (P97; see also P99). Most of the objects 
were followed photo-metrically for over a year. Except for a few SNe discovered before 1995, 
most spectra were obtained at an early enough epoch at the Keck 10-m telescopes to observe 
both the characteristic Type la SN spectrum and the host galaxy spectrum. Redshifts were 
determined primarily from narrow emission and absorption lines in the host galaxy spectra. 
In this paper we are mainly concerned with 35 of the 42 fully-analyzed SNe reported in P99 
which are listed in Table 1. The redshifts for this 35-supernova subset, 0.30 < z < 0.70, are 
such that "cross-filter K corrections" are made from the observed Kron-Cousins R to the 
rest-frame B band (Kim, Goobar, & Perlmutter 1996; Nugent, Kim & Perlmutter 2001), for 
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which lower-redshift data are available for comparison. Of the 42 SNe, 7 are omitted here 
as follows: Two z < 0.3 and three z > 0.7 supernovae, which have to be i^-corrected to the 
V and U bands respectively, are not included here. One distinctly faint outlier, SN19970, 
about 1.5 magnitudes fainter than the average, is also omitted. We have checked that it 
makes no significant difference whether this SN is included or not. SN1994G is omitted be- 
cause of sparse and late R photometry; it was observed primarily in the J-band. We consider 
one of the supernovae included in this study, SN1994H, to be unlikely to be Type la, based 
on our more recent analyses (P99, and Nugent, Kim & Perlmutter (2001)). It was included 
here for cross-comparison with the analyses of P97 and P99 (in which it was excluded from 
the primary analyses); we have checked that its light curve mimics a Type la SN quite well, 
and that including it in this study makes no significant difference. The lower-redshift sam- 
ple used here consists of 18 SNe selected from the 29 SNe of the Calan/Tololo set (Hamuy 
et al. 1995, 1996), based on the criterion that they were discovered prior to 5 days after 
maximum light. This homogeneous sample is chosen since it is the low redshift sample used 
in P99 which required SNe in the Hubble flow ( beyond z ~ 0.02) for comparison with the 
high redshift sample. We have also examined an additional set of 18 SNe which satisfy the 
criterion of discovery prior to 5 days after maximum, as discussed below. 

The SCP search strategy (P95) ensures early points by obtaining a set of several baseline 
"reference images" a few days after new moon, and another set of several "search images" 
about three weeks later. This ensures almost two weeks of dark time after discovery, for 
spectroscopic and photometric follow-up. Supernova candidates are found by differencing 
the pairs of images. The supernova may or may not be present, but in any case is much 
fainter on the reference set; this is established from archival references or references obtained 
a year later. The procedure yields one, or even two, image sets for most SNe on the early 
rising light curve. 



In the data analysis as described in P97 and P99, each SN's flux light-curve points are 
fitted to a i?-band template light curve using the nonlinear fitting program MINUIT (James 
& Roos 1994). For _R-band photometric data the function 



is fitted to the data by adjusting the intensity l max at maximum light, the time of maximum 
t max , the stretch factor s, so that the width factor w = s(l + z) scales the template time axis, 



3. Reduction to a Standard Z?-band Light Curve 




(1) 
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and a baseline level, b whose amplitude is found to be < 0.02 and allows for small corrections 
in the background galaxy subtraction. 

The function f R is generated from a £?-band template (in magnitude) (F B (t)) K- 
corrected to the i?-band for the given redshift, z, converted to flux and renormalized to 
unity at light maximum: /j?(0) = 1. We use as a template "SCP1997", the same used in P97 
and P99, and given in flux in Table 3 (Leibundgut 1988). SPC1997 was based on empirical 
fits by Leibundgut (shown in column 2 of Table 3 in flux ) , tabulated for —5 <t< 110 days 
Leibundgut (1988), and then extended by us with updated data to both late and early times. 
The constant slope in magnitude dmsldt = 0.0166/day for days 50-110 was extrapolated to 
even later times. A piecewise quadratic curve was adjusted to give a reasonable fit to the data 
for 8 nearby Type la SNe 16 for — 15 < t < — 5 days, and the slope defined by the two earliest 
points in magnitude (— 0.5041/day) was extrapolated to earlier times. This linear extrapo- 
lation in magnitude is exponential in flux. Only a handful of measurements were available 
for t < — 10 days, and the earliest behavior rested on single measurements of SN1994ae at 
— 12.9 days and SN1994D at —14.8 days in the observer system. This template was used for 
our cosmology analysis (P97, P99). Note that since all fits to data were performed starting 
with this template, this template has a stretch factor of s = 1 by definition. 

The results of the fits are shown in Table 1 which gives the redshift, z, the measured 
width factor, w, in the observer system and the corresponding stretch factor, s, as well as 
the the measurement errors, a w and a s , and the x 2 / DoF for the MINUIT fits. Table 2 gives 
the same information for the Calan/Tololo low-redshift data. 

The systematic uncertainty in the time-dependent K corrections, recalculated for P99, 
are estimated as < 0.02 magnitudes for the light-curve phases later than ~ —14 days. 
An effect we need to consider in determining the early time light-curve behavior is how 
uncertainties in the i^-corrections at these early epochs propagate into uncertainties on the 
fT-corrected 5-band fluxes. In particular, we need to examine all the points which occur 
prior to day —14 with respect to maximum light. Later than this epoch we have the necessary 
spectroscopy of a variety of SNe la to perform accurate if-corrections. Prior to this we have 
tied all the K-corrections to our knowledge of the spectral energy distribution (SED) of 
SNe la at this earliest epoch. 

At a redshift of z « 0.5, the uncertainty in the ^-corrections as a function of epoch is 
negligible due to the nice alignment of the B- and i?-band filters (Kim, Goobar, & Perlmutter 



SN1990N, SN1990af, SN1991T, SN1992A, SN1992bc, SN1992bo, SN1994D and SN1994ae 



6 



1996). For SNe la at redshifts below this value (to z = 0.3, our nominal cutoff for SN la 
included in this study) the B-to-R K-correction (Kbr) involves knowledge of SED redward 
of the 5-band filter. For SNe la beyond this redshift (to z = 0.7, our other cutoff) Kbr relies 
on our understanding of the [/-band behavior of SNe la. For a simple and very conservative 
test, we can gain an understanding of the uncertainties by looking at the i^-corrections 
at these two extreme redshifts, which involve the greatest extrapolation and hence would 
produce the largest error. 

For z = 0.3 and day — 14 we find Kbr = 0.51 for a SN la from a spectral analysis. 
This value corresponds to the value that a ~ 15,000 K blackbody would give at this redshift. 
From the models of Hoflich (1995), Nugent et al. (1997a) and Lentz et al. (2000), we can see 
that in the optical spectra redward of 4000A, at the earliest epochs after explosion, SNe la 
closely approximate blackbodies with a very minimal contribution from the lines (< 10% ) 
to the overall SED. The exact temperature is not critical in going from 15,000 K to 100,000 
K. Using a 100,000 K blackbody as a model of the earliest epoch for a SN la, (a gross 
overestimate of its temperature since the 7-rays from the decay of 56 Ni have not had enough 
time to propagate out and heat up the atmosphere) Kbr only changes by 0.07 magnitudes 
This is well below the the difference that can be distinguished by the data sets at these 
very early epochs. Only one SN la above a redshift of 0.6 (where the i^-corrections start 
to involve a decent amount of extrapolation), SN1995at, contributes to the early light curve 
prior to day —14. SN1995at was observed at day —16.3 and was at z — 0.655. The effect of 
this single photometry point on our calculations is negligible. In addition, all of the models 
mentioned above, show very little difference in the overall SED over the epochs —14 to —16.3 
days. 



3.1. A Composite Light-curve 

We now introduce the concept of a "composite light curve", constructed by linearly 
compressing or expanding the time axis for each SN such that all the low- z and high-z data 
points can be plotted on a single curve. To study the light-curve data in this composite 
form, a further step is added to the procedure of P97 and P99: After the time of maximum 
t max and the stretch factor for each SN is obtained from the fit, the individual i?-band data 
points are i^-corrected back to points in the equivalent rest-frame B band. 

Figure 1 presents all photometry points in composite light curves for the template 
"Parab-18" which is an improved template over SCP1997, as discussed below. The left-hand 
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figures (a), (c), and (e) show a data point for each night's observation of each supernova, 
whereas the right-hand figures (b), (d), and (f) show one-day averages over all supernovae. 
The K-corrected SCP points for SNe with 0.30 < z < 0.70 are shown as solid red circles. 
Shown as blue squares are -fT-corrected points [Kbb) for the 18 SNe from the Calan/Tololo 
study in the 5-band (Hamuy et al. (1996)) with z < 0.11. Figures 1 (a) and (b) show the 
£?-band photometry points in the observer system displaced to t — at light maximum, and 
normalized to unit intensity at t — 0. There are no corrections for stretch or width; times of 
observation relative to maximum light are used. In Figures 1(c) and (d) we have transformed 
the time axis from the observer frame to the rest frame by dividing by the appropriate 1 + z 
for each data point of each of the SNe. Finally, in Figures 1(e) and (f) for each supernova 
we also divide the time scale for each point by the fitted stretch factor s. By this stage 
essentially all of the dispersion has been removed, and the corrected points for both the SCP 
and Calan/Tololo data fall on a common curve at the level of the measurement uncertainty, of 
typically 2 — 4%. It is remarkable that application of this single stretch timescale parameter 
results in such a homogeneous composite curve. 



3.2. The Templates 

While we used the exponential rise in our fits to s, t max , I ma x and b with the SCP1997 
light-curve template, in reality there must be a definite explosion date t exp . To find an 
improved early template for the composite light curve we use a parabolic approximation to 
the initial light curve, and determine its two parameters, explosion date t exp and timescale 
parameter a, by calculating a least-squares fit of the function I(t)/I(0) = a(t — t exp ) 2 to the 
UN-averaged stretch-corrected composite flux data (which, unlike magnitudes, have Gaussian 
errors). (Since the fits were to the stretch-corrected composite data, the actual explosion 
day of any individual supernova is t = s x t exp .) We determined that fitting the early- 
time parabola to the later-time exponential template (matching both the value and first 
derivative) was not very sensitive to the join date in the epoch around -10 days. 

Using this additional form for the early part of the light curve, we can compare the 
combined composite light curve of Figures 1(f) with three forms of the template F B (t): 

1. The SCP1997 template: This light-curve template (listed in Table 3 ) begins with a 
linear rise in magnitude (exponential in flux); it was used in P97 and P99. The width 
factor and stretch factor values given in Tables 1 and 2 come from fits to this template 
and are supplementary to Tables 1 and 2 in P99 . 

2. The "Parab-18" template: This light-curve template (listed in Table 3 ) begins with a 
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parabolic rise in flux; its explosion date at —17.6 days is based on a parabolic fit to the 
composite data (starting from the SCP1997 template). A cubic spline fit to the data 
beyond the "join date" at —10 days improves the template fit as compared to SCP1997 
(or Leibundgut), in particular near -5 days. This template thus gives the best fit to the 
SCP data, although the fit to Parab-20 (see next item) is nearly as good. Whenever 
so stated the analyses in this paper are based on this template. 

3. The "Parab-20" template: This light-curve template (listed in Table 3) also begins 
with a parabolic rise in flux, but with an explosion date at —20.0 days based on the 
recent early data of Riess et al. (1999a). At later times, after day —10, it is the same as 
Parab-18, and thus differs by only a small amount from the SCP1997 template. This 
template also gives an acceptable fit to the composite data, as also found by Aldering, 
Knop, and Nugent (2000; hereafter AKN2000). 



Figure 2 a based on Parab-18, shows that this template accounts well for all the data 
points in our composite sample. Both the SCP and the Calan/Tololo data scatter about this 
template with a dispersion that is expected from their measurement uncertainties. From 
the residual plots of Figure 2 (b) and (d) it is clear that this day-by-day uncertainty is less 
than 2 — 4% (0.02-0.04 mag) for almost the entire period from day —10 to day 40. The \ 2 
distributions for the residuals is given in Table 4 (Parab-18) for the 35 SCP SNe in 5 day 
intervals from -25 to 40 days and, in the last line, for this entire interval. 

A good fit to the data is also obtained for Parab-20, corresponding to t exp = —20 days, 
for which we give the residual plot in Figure 2(d) and the x 2 distributions for the residuals 
in Table 5. Fitting the data to this new template, Parab-20, results in the fitting program 
MINUIT slightly re-adjusting stretch s, intensity at maximum light I ma x and the time of 
maximum t max for each SN. The average change in m e ^ ecttve from the values found using the 
SCP1997 template, for the 35 SCP SNe, is —0.005. The corresponding average change in 
m eJ.fective for thg lg Calan/Tololo SNe is -0.014. Thus the net change in difference between 
low and high redshift m e J fectwe is 0.009. This is the only quantity that enters into the 
cosmological measurements via the distance modulus. The rising part of the curve for both 
templates Parab-18 and Parab-20 fits our data equally well with only a negligible change 
in the overall x 2 /DoF from 0.905 to 0.906 for the epoch -25 to 40 days. Either of the 
above templates are an improvement over SCP1997 for which we obtain x 2 /DoF = 0.966 
(see Table 7). However within the accuracy of our data, any of the 3 templates is perfectly 
acceptable. 17 Note that the composite light curve depends in fine detail on the template 



17 Given the low-redshift data presented by Riess et al (1999a), we believe that the Parab-20 template is 
likely to be the best current estimate for a single common stretch-parameterized light curve. 



9 



used to fit and then match the peak and stretch of each individual SN's light curve. 

For the Parab-18 template we found a = 0.0085 ± 0.0010 and t exp = -17.6 ± 0.5 days; 
a preliminary value of —17.5 ± 0.4 was presented in the conference reports of Groom (1998) 
and Goldhaber (1998a, 1998b). However it is important to reckognize that these uncertainties 
reflected only the error in mapping the fixed composite light curve to a parabola after fitting 
with an assumed template (and thus did not include errors which depend on the shape of 
the early portion of the template used). Since the fits are sensitive to the early portion of 
the light curve, the total uncertainty on t exp is much larger. 

For comparison with AKN2000 we have analyzed the same more restrictive set of 30 
SCP SNe with 0.35 < z < 0.65 and leaving out SN1997aj. 18 A simple X 2 study of 
these SNe, fitting them to a series of parabolic-rise templates as a function of explosion 
dates, yields higher statistical uncertainties t exp = — (17.8^'g statistical) corresponding 
to a la range of —19.7 to —16.8 days. This statistical uncertainty is consistent with the 
more comprehensive analysis of AKN2000, which accounted for the full multidimensional 
probability distribution as well as upper limits on several systematic uncertainties, and found 
texp = — (18.3jl^;| statistical j^J systematic) at tj oin = —10 days. Note that the exact shape 
of the early portion of the light curve does not significantly affect the values or uncertainties 
of the cosmology results of P99 (see AKN2000 as well as our results above). 

Only future data will be able to tell whether all type la SNe can be adequately described 
by s only or whether some variation in t exp will also be needed. This will require very early SN 
detection, corresponding K-corrections, and spectra as is envisioned for the satellite project 
SNAP, currently under design. 19 

Table 6 gives the corresponding data for the 18 Calan/Tololo SNe from -10 to 40 days. 
Since there is no significant data before -10 days, Parab-18 and Parab-20 coincide in this 
case. For these data the first two epochs give a rather high x 2 / DoF of ~ 2. We note that 
at these two epochs the data come primarily from only two SNe; of these only SN1992bc 
has a poor overall \ 2 ( see Table 2) and in particular is responsible for the high \ 2 values at 
the two epochs. To improve the understanding of the early epochs we have also studied 18 
additional low redshift SNe which are however not discovered in a homogeneous search and 
are not all in the Hubble flow. These SNe consist of 16 of the 22 SNe given by Riess et al. 
(1999), which fulfill the criterion of being discovered prior to 5 days after maximum light, 



18 The fits for SN1997aj undergoe a large change in stretch, from 0.94 to 1.52, in going from the template 
SCP1997 to Parab-20. Because of this feature this SN was excluded by AKN2000 and in the 30 SN sample 
here However in the x 2 distributions quoted in Tables 4 and 5 all 35 SNe were included. 

19 See http://snap.lbl.gov for information on SNAP, the SuperNova Acceleration Probe. 
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as well as SN1990N from Lira et al. (1998) and SN1998bu from Suntzeff et al. (1999) as 
extended to earlier times (before —10 days ) by Riess et al (1999a). This low redshift sample 
fits the Parab-20 template from -15 to 40 days with a x 2 1 DoF of 1.015, thus agreeing well 
with the 18 Calan/Tololo SNe in the overlap region. This study of the second set of 18 SNe 
will be presented in a future publication. 

3.3. The Universality of the Stretch Factor 

The stretch factor applies both to the rising and to the falling part of the light curve. 
A qualitative impression of this can be obtained from a comparison of Figures 1(c) and (e) 
as well as (d) and (f). To obtain a quantitative result we took two approaches. In Table 8 
we show the x 2 distribution for the case in which only 1 + z but not s has been factored out 
of the observer frame distribution, this corresponds to Figures 1 (c) and (d). Here we note 
large x 2 values both before and after maximum light, indicating that ignoring s affects both 
distributions. Secondly we have taken the individual photometry points, of the composite 
curve and split them into two groups at epoch t = 0. We then fitted each half with MINUIT 
to determine a stretch. This gave s — 1.03 ± 0.04 with x 2 — 352 for 366 DoF for the epoch 
< t max and s = 1.00 ± 0.02 with x 2 = 894 for 996 DoF for the epoch > t max . The period 
> t max encompassed all the photometry points out to over 1 year. The equality of s, within 
statistical errors, for the two groups shows that the same value of s applies to both. 

As can be noted from Tables 4 and 6 as well as Figure 2 (b) and (c) the factor s combined 
with the factor 1 + z brings all the 35 SCP SNe and the 18 Calan/Tololo SNe in agreement 
with a single curve from -25 to over 25 days as well as can possibly be determined with the 
present datasets. The stretch factor parameterization appears to work quite well at least out 
to day 40, although the high-redshift data is less constraining at these late dates, and there 
is evidence of at least two low-redshift supernova that may vary from the template at these 
later epochs (see AKN2000). 

4. Test of Cosmological Light Curve Broadening and Check for SN Evolution 

These data also produce compelling evidence that the observed explosion of the su- 
pernova itself is slowed 20 by the factor 1 + z. This provides independent evidence for 



20 This light curve broadening has sometimes been described as "time dilation" , at variance with usage in 
physics texts (e.g. Weinberg (1972)). 
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cosmological expansion as the explanation of redshifts. Although this hypothesis has proven 
to be consistent with observation for over half a century, persistent doubts are still occasion- 
ally expressed (e.g. Marie et al. (1977); Chow (1977); La Violette (1986); Arp (1987); Arp 
et al. (1990, 1994); Narlikar & Arp (1993)). Surprisingly, until recently very few direct tests 
of this expansion have been performed. A test by Sandage & Perelmutter (1991), who cal- 
culated the surface brightness of brightest cluster galaxies over a range of redshifts, showed 
compelling evidence for expansion but did not reach a definitive conclusion due to possible 
systematic errors. An argument has also been made that some gamma ray bursts (GRBs) 
should be at cosmological distances because of the observation that for the "longer GRBs" 
the length of the bursts was inversely correlated to the brightness of the GRBs (Piran 1992; 
Norris et al. 1995). The discovery of GRBs at cosmological distances strengthens this argu- 
ment (see for example Metzger et al. (1997)); however since the intrinsic length of a given 
GRB is unknown, this remains a qualitative argument (Lee et al. 2000). 

Using supernova light curves to test the cosmological expansion was first suggested by 
Wilson (1939) (see also Rust (1974)). Over the last decade it has become clear that Type la 
SNe, found nearby and at cosmological distances, provide superb and precise clocks for such 
tests. We presented the first clear observation of the 1 + z light curve broadening, based on 
our first seven high z Type la SNe (Goldhaber et al. 1995). Leibundgut et al. (1996) later 
also gave evidence for the effect, using a single high-z supernova. More recently, Riess et al. 
(1997) showed evidence that the spectral features of Type la SNe can be timed sufficiently 
well to measure the time interval between two spectra taken 10 days apart in the observer 
system. Applying this method to one supernova gave results consistent with 1 + z light-curve 
broadening at the 96.4% confidence level. With the current dataset, we can now demonstrate 
the light-curve broadening with a larger, statistically significant sample. 

Light curve width factor w and stretch factor s versus 1 + z are shown in Figure 3(a) 
and Figure 3(b) respectively for the Calan/Tololo and SCP SNe. We are interested in the 
possible variation of w and s with z. We test for the z dependence of these distributions by 
fitting straight lines to the data. If the SNe in the distant past were different, i.e. due to 
evolution effects, the distribution might show a slope ds/dz. The fits described below were 
carried out for the entire 18 Calan/Tololo SNe together with 40 of the SCP SNe. Two SCP 
SNe, SN1992H and SN1992br, have stretch factors outside the range 0.7 < s < 1.3, and 
are excluded from the fits as outliers representing non-Gaussian tails to the distributions or 
aberrant objects. As discussed in P99, these will not be important for the present analysis. 

In performing the fit of the function s = a + bz to the data shown in Figures 3(b) the 
total uncertainty for each point is ay = a/°"s + $s where 8 2 s is the intrinsic stretch dispersion 
and a s from Tables 1 and 2 is from the individual light-curve fit uncertainty. We estimate 
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this intrinsic dispersion by requiring that the reduced x 2 is near unity. This yields an rms 
deviation S s ~ 0.08. The intrinsic width 5 S of the stretch factor distribution can also be 
obtained from a Gaussian fit to Figure 4 in P99. This gives a consistent value for 5 S of ~ 0.1. 
The errors used in the fit to the w distribution shown in Figures 3(a) are treated similarly. 
In Figure 3 (a) and (b) both errors a s and ay are shown on the error bars. 

We find that our data is consistent with 1 + z lightcurve broadening and a redshift 
independent SN stretch distribution. Fits to Figure 3 (a) and (b) yield different information. 
The linear fit shown in Figure 3(a) has dw/dz = 1.07 ± 0.06. If there were no z dependence 
in w this slope would be 0, here assuming no evolutionary change in s. Hence the evidence 
for the presence of a 1 + z factor is 1.07/0.06 ~ 18 standard deviations. Figure 3(b) shows 
a slope ds/dz = 0.05 ± 0.05. The extent by which this slope differs from measures the 
possible evolution effects on s, here assuming a 1 + z dependence of w. The result indicates 
that ds/dz < 0.09 at the 95% confidence level out to z ~ 0.8. We obtain essentially the 
same results from fits in which the 2 outliers are included (for a total of 60 SNe), and when 
the 35 rather than the 42 SCP SNe, less 2 outliers, are used (for a total of 51 SNe). Note 
that this analysis does not account for sample selection effects — such as possible preferred 
discovery of high stretch SNe near the flux limit of each SNe search — which may need to 
be accounted for in future datasets. 

When we compare this result to the alternate theories, it is clear that they are severely 
challenged or simply ruled out. The tired light theories (Zwicky (1929); Hubble & Tolman 
(1935); Hubble (1936); Marie et al. (1977); Chow (1977); La Violette (1986) would not yield 
this slowing of the light curves, and thus do not fit this dataset. Variable mass theories 
(Narlikar (1977); Narlikar & Arp (1993, 1997) would apparently require a series of coinci- 
dences: the radioactive decay times, the timescales for the radiative transport processes in 
the SN atmosphere, and the atomic transitions would all have to vary as (1 + z) due to 
variations in the masses of elementary particles for the theory to account for the data. All 
these effects would have to result in making the observed spectra similar. 

5. Discussion 

We have shown that a single stretch factor varying the timescale of the SNe la accounts 
very well for the restframe B band light curve both before and after maximum light (up 
to 40 days past maximum). It is not understood from the current state of the theory of 
SNe la whether this is a fortuitous coincidence or a reflection of some physics timescale 
(Hoflich 1995; Hoflich & Khokhlov 1996; Nugent 1997; Pinto & Eastman 2000; Arnett 2000; 
Hillebrandt & Niemeyer 2000). Although we do not necessarily expect this single timescale 
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stretch factor to hold for all wavelength bands, we do have evidence that it applies in the V 
band and (with less confidence) the U band. We will analyze these results in future work. 

The current analysis does not address the relationship between the SN la light curve 
stretch and its peak luminosity. This requires another analysis that will be presented else- 
where. However, it is important to note that the current B band template stretched by a 
factor s fits the data as well as any parameterization can, given the current data sets. It will 
therefore yield a peak-luminosity correlation with as small a dispersion as can be obtained 
by any other B-band light-curve parameterization. 

The comparison of the low-redshift and high-redshift composite supernovae light curves 
in Figure 2 provide a first-order test for evolution of supernovae as we go back to the z ~ 0.5 
epoch. The close match of the light curves is suggestive that little evolution has occurred. 

A recent analysis of low-redshift early light curves by Riess et al (1999a, 1999b) suggested 
that the explosion date of the low-redshift SNe la was earlier than that of the high-redshift 
SNe. The analyses of our early light curve uncertainty in Section 6 and in AKN2000 show that 
the difference between the two data sets is not very significant at this stage, less than ~2<r. 
This will be an interesting region of the light curve to pursue with future data, particularly 
over a range of host galaxy environments that would be expected to show variations in 
metallicity and progenitor ages. We are currently working with the supernova research 
community on generating a much more extensive low-redshift dataset to study this question. 
However, it is important to remember that - as pointed out by AKN2000 - the template 
differences have a negligible influence on the corrected peak magnitudes of the P99 SNe, and 
thus the cosmological parameters derived therein are unchanged. 

Finally, it is interesting to note that while the redshift of the light measures the expansion 
of the universe with a "microscopic clock" of period, typically T = 2 x 10~ 15 seconds, our 
"macroscopic clocks", the Type la SNe, measure the expansion over a ~ 4 week period or 
T 2.4 x 10 6 seconds. The 1 + z expansion effect is thus consistent for two time periods 
which differ by 21 orders of magnitude. 
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Table 1. Fit parameters for 42 Type 1a Supernovae reported 
by the Supernova Cosmology Project (Perlmutter et al. (1999)). 
For the present analysis we require good i?-band photometry and K 
corrections to B, so the 7 objects with qualifying comments are not 
used. 

Comments 
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* No -R-band data before 14 days; primary photometry from / band. 
** Here the V template was used while a template between V and R would be more 
appropriate. A fit to Parab-18 gives x 2 /DoF = 2.8. 
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Table 2. Fit parameters for the Calan/Tololo data 
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Table 3. B templates, in Normalized Flux, used by the Supernova Cosmology Project 

The headers: "SCP1997" refers to Leibundgut (1988) from day -5 on up, with a linear extension (in magni- 
tude) to earlier times shown in flux here. As discussed in the text, this is the template used in Perlmutter 
ct al. (1997, 1999). Tables 1, 2 and 7 as well as Figure 2(e), are based on this template. "Parab - 18" uses a 
parabolic turn-on at —17.6 days from a fit to the 35 SCP SNe. Tables 4 and 6 and Figures 1 and 2, except 
2(c) and 2(e), are based on this template. "Parab -20" uses the parabolic turn-on at -20 days following the 
early rise-time data of Riess et al. (1999a). Table 5 and Figure 2(c) are based on this template. After -10 
days the two parabolic templates coincide. 
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49 


0.055 


0.055 


12 


0.498 




0.511 


50 


0.054 


0.054 
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Table 4. y 2 Fits of the SCP Data to the 
Composite Parabolic Template "Parab-18" 
Corresponding to t exp = -17.6 Days 
Columns 1 and 2 define the Epoch Interval, 
Columns 3 and 4 give the x 2 an d 
X 2 per Degree of Freedom. See Fig 2 (b) 



Epoch 


x 2 


Y 2 /DoF 


Start 


End 






-25 


-20 


27.8 


1.16 


-20 


-15 


23.7 


0.99 


-15 


-10 


35.8 


0.58 


-10 


-5 


72.4 


0.95 


-5 





101.7 


1.02 





5 


111.6 


0.88 


5 


10 


74.5 


0.89 


10 


15 


62.2 


0.91 


15 


20 


25.6 


0.47 


20 


25 


83.3 


0.95 


25 


30 


70.1 


1.02 


30 


35 


49.1 


0.88 


35 


40 


25.8 


1.51 


-25 


40 


763.5 


0.905 



Table 5. x 2 Fits of the SCP Data to the 
Composite Parabolic Template "Parab-20" 

Corresponding to t exp = -20 Days 
Column headings as in Table 4. See Figure 2(d) 



Epoch x 2 X 2 /DoF 

Start End 



-25 


-20 


29.3 


1.01 


-20 


-15 


27.5 


0.65 


-15 


-10 


28.4 


0.71 


-10 


-5 


71.2 


0.94 


-5 





105.0 


1.00 





5 


115.6 


0.92 


5 


10 


70.1 


0.86 


10 


15 


63.6 


0.96 


15 


20 


26.1 


0.49 


20 


25 


84.6 


0.95 


25 


30 


67.9 


1.03 


30 


35 


52.0 


0.91 


35 


40 


21.4 


1.34 



-25 40 762.6 0.906 
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Table 6. \ 2 FlTS OF THE Calan/Tololo data 
to the Composite Parabolic Template* 



Column headings as in Table 4. See Figure 2(c) 



Epoch 


x 2 


X 2 /DoF 


Start 


End 






-10 


-5 


25.3 


1.95 


-5 





39.7 


2.09 





5 


17.1 


0.47 


5 


10 


59.5 


1.06 


10 


15 


55.7 


1.39 


15 


20 


28.9 


1.25 


20 


25 


14.8 


0.93 


25 


30 


11.2 


0.62 


30 


35 


21.1 


1.32 


35 


40 


17.0 


0.81 


-10 


40 


290.3 


1.143 



In this Table some of the photometry 
errors were extremely small. We used 

a minimum error of 0.007, in 
normalized units, here. 



Table 7. \ 2 Fits of the SCP Data to the 
Composite Exponential Template "SCP1997" 
the Extended Leibundgut Template 
Column headings as in Table 4. See Figure 2(e) 



Epoch 


x 2 


X 2 /DoF 


Start 


End 






-25 


-20 


27.5 


1.15 


-20 


-15 


24.3 


1.01 


-15 


-10 


36.8 


0.59 


-10 


-5 


105.2 


1.38 


-5 





105.7 


1.06 





5 


110.8 


0.87 


5 


10 


75.3 


0.90 


10 


15 


73.4 


1.08 


15 


20 


29.3 


0.54 


20 


25 


83.8 


0.95 


25 


30 


69.1 


1.00 


30 


35 


49.2 


0.88 


35 


40 


25.4 


1.49 


-25 


40 


815.7 


0.966 
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Table 8. \ 2 Fits of the SCP Data to the 

Parabolic Template "Parab-18" 
with only 1 + z but not s factored out 

This Table illustrates the effect of the 
factor s both before and after Maximum light. 
Column headings as in Table 4. 



Epoch x 2 x 2 /DoF 

Start End 



-25 


-20 


24.2 


1.86 


-20 


-15 


30.3 


1.08 


-15 


-10 


231.8 


4.55 


-10 


-5 


250.2 


3.05 


-5 





104.3 


1.03 





5 


112.4 


0.89 


5 


10 


125.6 


1.38 


10 


15 


57.9 


0.95 


15 


20 


110.1 


1.84 


20 


25 


98.6 


1.15 


25 


30 


105.0 


1.36 


30 


35 


45.9 


1.12 


35 


40 


35.4 


1.18 



-25 40 1330.7 1.582 
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-40 -20 20 40 -20 20 40 60 

Observer system days from peak Observer system days from peak 




Rest system days from peak Rest system days from peak 



Fig. 1. — (a) The photometry points for the 35 SCP (full red circles) and 18 Calan/Tololo 
SNe (blue squares), fitted to Parab-18 with the maximum flux normalized to unity and the 
time of maximum adjusted to zero in the observer system, (b) shows the same data as in 
(a) averaged over one-day intervals and over each set of SNe. (c) and (d) show the same 
data, transformed to the rest system. In (e) and (f) the time axis for each photometry point 
is additionally divided by the corresponding stretch factor s. 
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Fig. 2. — (a) The same distribution as in Figure 1(f) with the SNe in the two data sets fitted 
to Parab-18. The template curve for Parab-18 is shown in black, (b) Residuals for the SCP 
(solid red circles) photometry fits to Parab-18. Here the epoch between -40 and -20 days 
has no SN signal, the data comes from the reference images which have shorter exposures 
and hence larger errors (see Table 4). (c) Residuals for the Calan/Tololo (blue squares) 
photometry fits to Parab-18 (see Table 6 ). (d) The open red circles show the residuals for 
the SCP data from the fit to Parab-20 (template not shown), where the early rise is based 
on the data of Riess et al. (1999a). See Table 5. (e) Residuals for the SCP data (open red 
diamonds) fitted to the template SCP1997 (template not shown). Note that the "bumps" 
near -5 days and 15 days (see also Table 7) were eliminated for the templates Parab-18 and 
Parab-20. 
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Fig. 3. — (a) The observed light-curve width factor w vs 1 + z. The blue squares correspond 
to the Calan/Tololo SNe, the solid red circles are for the subset of 35 SCP SNe used in this 
paper. The open red circles are for the remainder of the 42 SCP fully-analyzed SNe. The 
band delineated by the black dash dotted lines corresponds to stretch values 0.7 to 1.3 which 
encompass the bulk of the data, except for two outliers. The green line is the best linear fit 
to the data as discussed in the text. The band delineated by the two green dashed curves 
correspond to the ±lcr values, (b) The stretch s vs 1 + z. Stretch is defined as the observed 
light-curve width w divided by 1 + z for each SN. The points and lines are defined as in (a). 



